Comments on the entropy of seismic electric signals under time reversal 
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We present recent data of electric signals detected at the Earth's surface, which confirm the earlier 
finding [Phys. Rev. E 73, 031114 (2006)] that the value of the entropy in natural time as well as its 
value under time reversal are smaller than that of the entropy of a "uniform" distribution. Further- 
more, we show that the scale dependence of the fluctuations of the natural time itself under time 
reversal provides a useful tool for the discrimination of seismic electric signals (critical dynamics) 
from noises emitted from manmade sources as well as for the determination of the scaling exponent. 
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In a time series comprising N events, the natural time 
Xk = k/N serves as an indexfj [3] for the occurrence 
of the fc-th event. In natural time analysis, the time 
evolution of the pair of the two quantities (XkiQk) is 
considered, where Qk denotes in general a quantity pro- 
portional to the energy released during the fc-th event. In 
the case of dichotomous electric signals (SES) activities 
(i.e., low frequency < 1Hz electric signals that precede 
earthquakes, e.g., see Refs.[3j, U, |5j, |6J]) Qk stands for the 
duration of the fc-th pulse. The entropy S in natural time 
is defined^ as 



S = (%lnx) - (x) ln(x) 



(1) 



where (/(%)) = Y,k=iPkf{Xk) and p k = 
The value of the entropy upon considering the time re- 
versal T, i.e., Tpk — PN-k+i, is labelled by S-. 

SES activities (critical dynamics) exhibit infinitely 
ranged long r ang e temporal correlations [1, 0, [1(3] which 
are destroyed |1C| after shuffling the durations Qk ran- 
domly. An interesting property emerged from the data 
analysis of several SES activities refers to the fact@ 
that both S and S— values are smaller than the value 
of S u (= 1/2 In 2 - 1/4 « 0.0966) of a "uniform" distribu- 
tion (defined in Refs. pi @, EH, EH , e.g. when all pk are 
equal), i.e., 



S — *\ Six 



(2) 



This finding, which does not hold [13] for "artificial" 
noises (AN) (i.e., electric signals emitted from manmade 
sources), has been also supported by numerical simula- 
tions, e.g., in fractional Brownian motion (fBm) time 
series Q that have an exponent a. r> fa resulted from the 
Detrended Fluctuation Analysis (DFA)[l3, 0]| close to 
unity. We clarify that fBm (with a self-similarity in- 
dex H ss 1) has been found [161] as an appropriate type 
of modelling process for the SES activities. Thus, it 
seems [9j that the validity of the relation (2) stems from 
infinitely ranged long range correlations. The scope of 
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this paper is twofold: First, to provide the most recent 
experimental data that strengthen the validity of the re- 
lation (2), and second to point out the usefulness of the 
study of the fluctuations of the natural time itself under 
time reversal. 

Figure 1 depicts an electric signal, consisting of a num- 
ber of pulses, that has been recorded on November 14, 
2006 at a station labelled)!?]] PIR lying in western Greece 
(close to Pirgos city). This signal has been clearly col- 
lected at eleven measuring electric dipoles with electrodes 
installed at sites that are depicted in a map given in 
Ref.[l7j]. The signal is presented (continuous line in red) 
in Fig. 1(a) in normalized units, i.e., by subtracting 
the mean value and dividing by the standard deviation. 
(As for the actual amplitude -in mV/km- of this SES 
activity [HI, it is comparable to the one observed at the 
same station before the magnitude M rj 7.0 earthquake 
that occurred on Jan 8, 2006, see Ref.[l9(.) For the 
reader's convenience the corresponding dichotomous rep- 
resentation is also drawn in Fig. 1(a) with a dotted (blue) 
line, while in Fig. 1(c) we show (in red crosses) how the 
signal is read in natural time. The computation of 5* and 
S- leads to the following values: S = 0.070 ± 0.012, 
S- = 0.051 ± 0.010. As for the variancep], i, 0, 
Ki = (x 2 )~(x) 2 , the resulting value is m = 0.062±0.010. 
These values more or less obey the conditions K\ « 0.070 
and S, S- < S u that have been found to hold for other 
SES activities @|. 

A closer inspection of Fig. 1(a) also reveals the fol- 
lowing experimental fact: An additional electric signal 
has been also detected (in the gray shaded area of Fig. 
1(a)), which consists of pulses with markedly smaller am- 
plitude than those of the SES activity discussed in the 
previous paragraph. This is reproduced (continuous line 
in red) in Fig. 1(b) in an expanded time scale and for the 
sake of the reader's convenience its dichotomous repre- 
sentation is also marked by the dotted (blue) line, which 
leads to the natural time representation shown (dotted 
blue) in Fig. 1(c). The computation of S and 5_ gives 
S = 0.077±0.004, S*_ = 0.082±0.004, while m is found to 
be ki = 0.076 ± 0.005. Hence, these values also obey the 
aforementioned conditions (k% « 0.070 and S, S- < S u ) 
for its classification as an SES activity. 
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FIG. 1: (color online) (a) The electric signal recorded on 
November 14, 2006 at PIR station (sampling rate fex P ~ 
1Hz). The signal is presented here in normalized units (see 
the text). The corresponding dichotomous representation is 
shown with the dotted (blue) line. The gray shaded area 
shows an additional signal (consisting of pulses of smaller du- 
ration) superimposed on the previously mentioned signal, (b) 
The signal belonging to the gray shaded area in (a) is given 
here in an expanded time scale, while its dichotomous repre- 
sentation is marked by the dotted (blue) line, (c) How the 
signals in (a) and (b) are read in natural time: the continuous 
(red crosses) and dotted (blue asterisks) bars correspond to 
the durations of the dichotomous representations marked in 
(a) and (b), respectively. 
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FIG. 2: (color online) The log-log plot of the fluctuations A\i 
of the natural time under time reversal versus the scale I for 
fBm (a) and fGn (b). Time series of length 2 x 10 4 have been 
used. 
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FIG. 3: (color online) The log-log plot of A\i versus the scale I 
for three SES activities (Kl, K2 and U) and seven AN (nl-n5, 
n7 and nlO) treated in Ref. [13| (cf., these signals have enough 
number of pulses in order to apply the present analysis). The 
values of A\i are divided by the corresponding values A^s at 
the scale I — 5. 
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We now proceed to the second goal of this paper. The 
way through which the entropy in natural time captures 
the influence of the effect of a small linear trend, has 
been studied in Ref.0]. Namely, the "uniform" distri- 
bution, p(x) — 1, where p{x) is a continuous probabil- 
ity density function (PDF) corresponding to the point 
probabilities pk used so far, has been perturbed by a lin- 
ear trend e in order to obtain the parametric family of 
PDFs: p(x; e) = l + e(x- 1/2). Such a family of PDFs 
shares the interesting property Tp(x; e) = p{x\ ~ e ): i- c i 
the action of time reversal is obtained by simply chang- 
ing the sign of e. It has been shownQ that the entropy 
5(e) = 5[p(x;e)], as we U as that of the entropy under- 
time reversal S-(e) = S[Tp{ X ;e)], S-(e) = S(-e), de- 
pend non-linearly on the trend parameter e: 



If we consider the fluctuations of this simple measure 
upon time-reversal, we can obtain information on the 
long-range dependence of Qk- Actually, we shall show 
that a measure of the long-range dependence emerges 
in natural time if we study the dependence of its vari- 
ance under time-reversal Ax 2 = E[(x — Tx) 2 ] on the 



window length I that is used for the calculation. 
Tp k = pi-k+i, we have 



Since 



S(e) 



e 

72 



-W 1 

12 J \2 
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(3) 



However, it would be extremely useful to obtain a linear 
measure of e in natural time. Actually, this is simply the 
average of the natural time itself: 



(x) = 



(4) 



A X 2 = E[( X - T X ) 2 } = E ■ 



1 k 

,fe=i 



■pi-k+i) 



( 5 ) 

where the symbol E[. . .] denotes the expectation value 
obtained when a window of length I is sliding through 
the time series Qk- By expanding the square in the last 
part of Eq. ([5|) , we obtain 



Ax 2 = (y ) -Pl-k+l) 2 } + Yl -fT E l(Pk ~ Pl~k+l){p m -Pl- m +l)]. (6) 

k=l ^ ' k=£m 

The basic relation that interrelates pt isfTT| 

i 

fc=i 

or equivalently pk = 1 — Em#P™- By subtracting from the last expression, its value for k = I — k + 1, we obtain 
Pk - pi-k+i = -J2 m =£k(P™ ~ Pl-m+l), and thus 

(p k -Pl-k+l) 2 = - E (P k - Pl-k+l)(Pm - Pl-m+l)- (8) 

m^k 

By substituting Eq.® into Eq.©, we obtain 

7 ) X/ E ^ Pk ~ Pl-k+l)(Pm ~ Pl-m+l)] + E "72~ E [(P fc -Pl-k+l)(Pm ~ Pl-m+l)] (9) 
n- — x m^k k^m 

i 



which simplifies to 

Ax 2 = - E ^ E t( Pfc -Pl-k+l)(Pm - Pl- m +l)] 

k,m 

(10) 

The negative sign appears because (pk — Pi-k+i) and 
(Pm — Pi-m+i) are in general anti-correlated due to 



Eq.©. Equation (JTUJ) implies that Ax 2 measures 
the long-range correlations in Qk'- If we assume that 
-E[(p fc -p ; _ fc+ i)(p m —Pi-m+i)] {k-m) 2xH /I 2 (cf. p k 
scales as 1/7, e.g. see [HI), we have that Ax 2 oc 1 4+2xh /I 4 
so that Ax oc l XH (where xh is a scaling exponent). 
In order to examine the validity of this result in the 
case when Qk are coming from a fractional Gaussian 
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noise (fGn) or fBm, we employed the following proce- 
dure: First, we generated fBm (or fGn) time-series X/. 
for a given value of H using the Mandelbrot- Weierstrass 
function (20l. |2ll |22| as described in Ref.0]. Second, since 
Qk should be positive, we normalized the resulting Xk 
time-series to zero mean and unit standard deviation and 
then added to the normalized time-series Nk a constant 
factor c to ensure the positivity of Qk — Nk + c (for the 
purpose of the present study we used c = 10.0). The 
resulting Qk time-series were then analyzed and the fluc- 
tuations of Axi versus the scale I are shown in Figure 2. 
We observe that for fBm the estimator \h ~ H, whereas 
for fGn xh ~ H — 1. The physical meaning of the present 



analysis was further investigated by performing the same 
procedure in the time-series of the durations of the elec- 
tric signals analyzed in Ref.[l3[. The relevant results are 
shown in Figure 3. Their inspection interestingly indi- 
cate that all seven AN correspond to descending Axi 
curves versus the scale I, while the three SES to ascend- 
ing curves. This fact is in essential agreement with the 
results obtained in Ref. 0] , which showed that the SES ac- 
tivities have a stronger memory than AN. In other words, 
apart from being an estimator of the scaling behavior, the 
ascending or descending behavior of Axi versus the scale 
I seems to provide a useful new tool for the classification 
of a signal as SES activity or AN, respectively. 
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